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There has been considerable debate on the onset of exotic spin phenomena in quantum wires due 
to enhanced many-body effects caused by the one-dimensional (ID) alignment of charge carriers. 

We explain various observed spin effects, such as a carrier density-dependent spin-flip in dilute 
quasi-lD systems and the variability of the spin polarization in quantum point contacts, by using 
an unrestricted Hartree-Fock approach with a three-dimensional (3D) Coulomb interaction. The 
model dimensionality is critical in identifying a complex pattern of magnetic phases varying with 
magnetic field and confinement. In the limit of vanishing magnetic fields, we show the emergence 
of a degenerate excited state with opposite spin polarization above a confinement-dependent ID 
concentration threshold, which is consistent with observations of a conductance plateau at half the 
conductance quantum Go/2 = e^//i, even in the absence of spin-orbit interactions. 


Quantum wires (QWRs) are nanostructures characterized by two-dimensional (2D) confinement exhibiting electronic 
modes transverse to the one-dimensional (ID) motion of charge carriers. The existence of these transverse modes 
has profound consequences on the interaction between carriers and crystal dynamics, as well as amongst carriers 
themselves, which uncovers a flurry of exciting properties. In the former case, carrier scattering undergoes size 
effectsi that affect the transport properties with important technological consequences for device electronics^. In the 
latter case, the ID alignment of interacting particles enhances Hartree and exchange interactions giving rise to exotic 
phenomena such as the formation of Wigner localization^ii^ or Luttinger liquids^i^. During the last two decades, a 
wide range of experiments stimulated by the observation of a transport anomaly in semiconductor quantum point 
contacts (QPCs)i“— and the demonstration of separate spin and charge excitations in have suggested the 

existence of spin-related transport or spontaneous spin polarization in ID system s as opposed to a Kondo-like 
effect due to the presence of a quasi-bound state^ '^^’^^ . A spin-polarized ground state, while forbidden in strictly 
ID systems by the Lieb-Mattis theoreruiSi^ could be achieved in realistic QWRs and QPCs, since these devices have 
two- and three-dimensional structures, potentially giving rise to additional phenomena. This has spurred interest in 
fully-electric manipulation of spin properties in quantum wires and point contacts^^^. 

In this article, we show that ID systems in longitudinal magnetic fields and in the absence of spin-orbit interaction 
(as, for instance, in a GaAs wire) can sustain a hierarchy of spin configurations depending on carrier concentration, 
energy and confinement. Specifically, as carrier concentration increases, the electron system in its ground states 
evolves from a fully spin-polarized state to an unpolarized state with a spin flip. In the limit of vanishing magnetic 
fields, there exists a concentration threshold above which the electron system exhibits an excited state with degenerate 
opposite-spin polarization. We also show that this threshold varies as the strength of the confinement is changed. The 
latter feature is particularly important for the technological application of electrostatically-confined wires for which 
the spin polarization can be controlled by electrical gating. 

Let us consider a wire with its axis of symmetry along the cc-axis. Electrons are confined along the y- and z-directions 
(perpendicular to the wire) by means of a potential Dconf (y, z), which we model as a superposition of two parabolic 
wells, i.e. Uconf {y,z) = ^m*ujyy^ Here, Uy and u)z are the strengths of the confinement along the y- 

and z-directions, respectively, and m* is the electron effective mass. The wire is placed in a magnetic field B = Bqx 
parallel to the axis of the wire, with an associated Zeeman term Uz = gUsBxO'- (<7 is the effective electron g-factor, 
PB — 2 ^ i® the Bohr magneton in the wire, q is the electron charge, and a is equal to -I- 1/2 or —1/2 for spin-up or 
spin-down, respectively.) For the sake of simplicity, the spin-orbit interaction is neglected. 

The use of the unrestricted Hartree-Fock model in the extreme quantum limit (i.e. when only one subband is 
populated)^ results in the expression for the energy E {kx , cr) of an electron in terms of its momentum kx and spin 
a (see Supplementary Methods), which reads: 
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Here, cob = is the cyclotron frequency, Uei is the Hartree term (which accounts for Coulomb repulsion amongst 
electrons and is proportional to the total concentration no) and [Axch is the exchange term. Ugi and C/exch depend 
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on the overlap function (^ab (p) which, in turn, is determined by the strength of the lateral confinement and the shape 
of the electron wavefunction. (See Supplementary Methods.) The pair of two integral equations for and 

E {kx,i) described by Eq. H]) is then solved to yield the spin-dependent concentrations Ua {n^, ni) for a fixed total 
electron concentration uq = n-f- -I- uj, . 

At zero temperature, both rif and nj, are associated with a (positive) spin-dependent Fermi wavevector = Trno-, 
so that the Fermi energy Ef = E {kf(^^)) = Ef { 71 ^). Then, after setting Ef (n^) = Ef (nj,), the relation between n-j- 
and ni can be written in terms of no and An^ = Ua — noj2, i.e. 
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This equation is discussed in more detail in the Supplementary Methods and holds only for 0 < |Ano-| < no/2, since 
0 < Ua- < riQ. If Eq. ([2]) yields the solution An.|-^) > no/2, the identity fc/( 4 ,(-|-)) = is no longer valid because 

n^(l) would then be negative, and thus must be set to zero in order to still satisfy Eq. ([j). In this case, nj^(.|-) = no, 
which corresponds to full spin polarization. 




Figure 1. Spin-polarized regimes, (a) Polarization P = (n-|- — nj,) /no as a function of the total concentration no at = 1 T 
and T = 0 K for non-interacting electrons (solid line) and the three possible polarization configurations for interacting electrons: 
d (dashed), f (dotted) and d* (dash-dotted). 


In order to illustrate our model, we consider a GaAs quantum wire with |(jf| = 0.44, m* = 0.067too and e = 12.9e(j2^, 
and where the transverse confinement strength is set to Hujy = hwz = 2.0 meV. Figure [T](a) shows the polarization 
P = (n-[ — ni) /riQ for interacting and non-interacting electrons, as a function of the total concentration ng, at = I T 
and zero temperature. The solid line in Figure[T](a) corresponds to the spin polarization for non-interacting electrons, 
i.e. Cab (p) = 0 in Eq. ([J). In this case P is negative because of the Zeeman interaction, Uz = gpsBxO', which lowers 
the potential energy of spin-down electrons, and hence nj, exceeds n-]- at the same Fermi energy. For concentrations 
n-o < ns = 0.817 x 10®cm“^, there is complete spin-down polarization and P = —1. For rig exceeding ub, one gets 
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(as described in the Supplementary Methods), so P = — and 


-1 < P < 0, 


which corresponds to partial spin-down polarization. As ng increases to infinity, P slowly approaches zero. This is due 

to the fact that the kinetic energy term 7/;(cr) = 2 m‘'’ ~ ’ "''^hich increases with rio-, contributes predominantly 

to the total energy at high concentration since Ez is independent of ng. As a result, the Zeeman splitting induced by 
Ez becomes insignificant. 

The dashed line in Figure (TJa) is one of the solutions to Eq. ([T]) for interacting electrons, which we call the “f” or 
“up” configuration, as P > 0. For this “t” configuration, the wire is fully spin-polarized when ng < = 

4.51 X 10^cm“^ and ng > ^ ^ io®cm“^. These limiting values are obtained by solving Eq. ([2]), 

setting An.|- = Ptiq/^. Between those two concentrations, there is partial spin polarization, with a minimum P = 0.32 
at ng = 0.54 x 10^cm“^. Spin polarization is opposite to the non-interacting case because of the existence of the 
exchange interaction, which lowers the zero-point energy of the ID energy subband. At very low ng, the exchange 
energy Pexch dominates the kinetic energy Tx, but since Ef (n^) = Ef {n^), rif increases and nj, decreases (otherwise, 
C4xch(4,) would be much more negative than Pexch(t)i which would cause Ef (nj,) to drop far below Ef (n-f-)). 
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For no > n™®®* (= 1.94 x 10^cm“^), two additional configurations emerge, both with P < 0: “4,” or “down” (dotted 
line in Figure [IJa)) and “J,*” or “down-star” (dash-dotted line). In both cases, P = —0.62 at the threshold concen¬ 
tration n™®®*, but the two configurations differ as no increases. In the “J,” regime, the polarization becomes stronger 
(i.e. more negative) until the wire is fully polarized (P = —1) above n-Q”" = 2.27 x 10®cm“^ (a value also obtained 
from Eq. @ with Anj, = -|-no/2). Meanwhile, in the “J,*” regime, spin polarization is weakened and tends to zero 
(approaching the Zeeman splitting for non-interacting electrons) as ng goes to infinity. Here we point out that the 
specihc numerical values of the concentrations Uq"®®*, etc., as well as the energies displayed on figures [T]131 are 
all confinement-dependent and will change for different confinement strengths (see Fig. 0. 

The different high-concentration behaviors of the three configurations (P = 1, P = — 1 and P —>■ 0 for the t, 4' and 4-* 
configurations, respectively) result from the interplay between the kinetic energy T^, the exchange interaction Pexch 
and the Zeeman splitting Uz to satisfy the condition Ef (ni) = Ef (nj,). At high no, Pexch dominates Uz, which can 
be attenuated in three ways: (I) If is smaller than n-[, which prevents the combined contribution Pexch( 4 .) + Uz{i) 
from becoming too negative; this leads to full spin-up polarization (f configuration). (2) If n^ is instead larger than 
n-f, which reduces the kinetic energy term T!c(t) relative to thus lowering Ef (n-|.) and resulting in full spin-down 

polarization (4, conhguration). (3) If nf is almost equal to n-f (i.e. P —>■ 0, just like for non-interacting electrons) in 
order to balance 74,(4,) s-iid since dominates Pexch(and Uz) at high concentrations; this is the behavior of the 
4,* configuration. 

Whether the electrons in the ground state of the system are in the f, 4' or I* configuration depends on which of 
these three configurations has the lowest energy. Figure [T](b) shows the average energy per electron as a function 
of no for all three possibilities. If ng is above tiq"®®*, when all three solutions are possible, the 4-* configuration has 
the lowest energy, with an energy difference with the other two solutions that grows as ng increases. While the 4- 
and 4.* configurations both emerge with the same energy Pavg = 5.2 meV at Uq"®®*, the separation between these 
two grows up to 0.7meV at no = 3 x 10^ cm“^. The f configuration has the highest energy, exceeding that of the I 
solution by a constant value of 0.4 meV. However, below the n ™®®‘ threshold, 4* is forbidden, so the system assumes 
the t configuration. As a consequence, at n™®®*, the ground state changes abruptly from a positive to a negative 
polarization. This spin-flip process is shown on the insets of figure HKb) . The polarization Pgnd of the ground state 
(top-left inset) drops suddenly from -|-I (t configuration, dashed line) to —0.6 (4.*, dash-dotted line) at n™®®*. The 
bottom-right inset displays the ratio najng for the ground state of the system. Below n™®®*, the system is in the 
f-configuration where n-f (solid line) exceeds (dashed line). In this region, there is full spin-up polarization for 
concentrations below 4.51 x I0^cm“^ and above 1.92 x I0^cm“^. Between these two values Uf reaches a minimum at 
no — 0.54 X 10® cm“^ where only two-thirds of the electrons in the wire have spin-up. At ng"®®*, the system switches 
to the 4 .* configuration, for which > n^. At this spin-reversal point, three-fourths of the electrons have spin down, 
but as no increases the difference between and n^ decreases. 



Figure 2. Magnetic field dependence of the spin-polarized regimes, (a) Magnetic phases diagram of the distinct 
spin-polarized regimes as a function of no and Bx'- (I) t configuration only, partially polarized; (II) t configuration only, fully 
polarized; (III) 4 (partial), 4 (partial), 4* (partial); (IV) 4 (full), 4 (partial), 4* (partial); (V) 4 (Ml), 4 (full), 4* (partial), (b) 
Polarization vs. no for different magnetic field strengths. In all cases, hwy = hu)z = 2.0 meV and T = OK. 


Figure [2Ka) displays the various spin polarization regimes in the wire as both the magnetic fields and the concentra¬ 
tions are varied. Only the 4 configuration is present in regions (I) and (H), with either partial or full polarization, 
respectively. These two regions are separated by a solid, almost parabolically-shaped line. For > 4.22 T the 4 con¬ 
figuration is completely polarized regardless of the concentration. The black dashed line corresponds to the minimum 
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concentration Uq"®®* for the emergence of the i and J,* configurations. In regions (III) and (IV) the f configuration 
is, respectively, partially and fully polarized. Finally, region (V) lies above the topmost solid line that indicates the 
minimum concentration for which the configuration is fully polarized. In this region only the 4-* configuration 
is partially polarized. 

The presence of different polarization regimes in Fig. [2]is consistent with direct measurements of the spin polarization 
in quantum point contacts^^ showing their variability as the magnetic field and concentration are changed. Further¬ 
more, the sudden polarization reversal predicted in our model (Fig. [Ub)) also explains the abrupt rearrangement of 
the spin-up and spin-down levels under a strong in-plane magnetic field as observed by Graham et al^^ and previously 
interpreted as an exchange-driven magnetic phase transition^^Ti^. Indeed, at Uq"®®*, the polarization of the ground 
state changes from positive to negative, so the energy of spin-up electrons suddenly rises above that of spin-down 
electrons, leading to the observed depopulation of the spin-up subband. 

Figure HKb) shows the effect of vanishing magnetic fields on the different spin polarizations in the wire, which evolve 
from a complex of configurations (f, 4^, and 4-*) for Bx ^ 0 to a situation that sustains an unpolarized state co-existing 
with two degenerate and symmetrically spin-polarized regimes at high concentrations for Bx = 0. One can clearly see 
that the t spin polarization at low concentration, and the 4-* spin polarization at high concentration collapse to an 
unpolarized configuration, while the 4- and t configurations are symmetric relative to each other. At the same time, 
the concentration threshold Uq"®®* decreases to a lower common value for up- and down-spin polarization. As will be 
seen in Fig. EKb), these two configurations are degenerate, characterized by a single Fermi level and equally probable, 
which is consistent with the absence of net magnetic moment in the wire. 



Figure 3. Spin-polarized regimes at Bx = 0. (a) Spin-polarized concentrations n-f and ni vs. no for the 4 configuration 
when Bx = 0. Inset: polarization of the 4- configuration vs. no. (b) Fermi energy vs. no for both the polarized (solid) and 
unpolarized (dotted) regimes. The vertical dash-dotted line marks the concentration for full polarization (no“" = 2.lx 10® cm“^). 
tkoy = hwz = 2.0meV and T = OK. 


The existence of spin polarization degeneracy when Bx = 0 can be shown from Eq. m where, asides from the trivial 
solution Aria = 0 that corresponds to a spin-unpolarized electron density, one obtains the approximate solution for a 
non-zero Aria in the limit |Ancr| <C ng (see Supplementary Methods): 


Aria 


±4 


C(7rno) . 




Ca b (tTUq) 
27r^ 


(3) 


Here, a* = 

’ q-^m* 


is the effective Bohr radius and Cab (’’’’^o) = 


d^Cabip) 


dp^ 


p—7Tno 


. This solution is only valid for no > 


’^o(Bx=o) ’ where Ug^'^®*^p^ satisfies the identity obtained by setting Aua = 0 in Eq. ([3|) : 
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Eigure[2Ka) is a plot of the spin-polarized concentrations in the wire for the 4- configuration {n\, > n^) when Bx = 0. 
The solid line corresponds to the unpolarized case {n^ = n-\- = no/2). The spin-polarized regime, represented by the 
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dashed curves, emerges at ng = n™®®* (= 1.54 x 10®cm“^). Above this threshold continues to increase with ng 
until the wire becomes fully polarized (nj, = ng, n-f = 0) for ng > (= 2.1 x 10®cm“^). The inset shows the 

polarization corresponding to this configuration as it changes from 0 at uq = Uq"®®* to — 1 at ng = 

Both spin-polarized configurations (f and J,) have the same energy, which is generally higher than that of the unpo¬ 
larized regime. This is shown in Fig. [3Jb), which is a plot of the Fermi energy as a function of total concentration. At 
^ons®t Fermi energy is 6.97meV for all configurations, but as ng increases the energy of the polarized configura¬ 
tions rises more rapidly, so at ng = the difference in energies between the two cases is about 1.5 meV. Therefore, 
at T = OK, the ground state of the electrons in the wire will be unpolarized, in agreement with the Lieb-Mattis 
theorem^!. However, at higher temperatures the system can be excited to one of the spin-polarized regimes (which 
are indistinguishable from each other in the absence of a H-field), especially for concentrations close to n™®®‘ when 
the energy difference between the polarized and unpolarized configurations is lower. 





Figure 4. Dispersion relation. Energy vs. wavevector when = 0, for three different concentrations: (a) just below 
the threshold for polarization; (b) just above the threshold; (c) significantly above the threshold, hjjy = hwz = 2.0 meV and 
T = OK. Ao (A:), Aj-jg) and A:/(o) are, respectively, the energy and the Fermi wavevector in the unpolarized configuration, while 
Ea (k) and fcj-(£,)are the energy and the Fermi wavevector in the polarized configuration for electrons of spin a (=t, 4)- 


The energy-momentum relation for spin-up and spin-down electrons at = 0 is illustrated in Figures|3Ka-c) for three 
different concentrations in the ], configuration (P < 0). In all cases there is an inflection point at the Fermi wave vector 
caused by the exchange interaction, which indicates the presence of a maximum in the carrier velocity, a minimum 
in the ID density of states, and an infinite effective mass. In Fig. |D(a), for which ng < n™®®* (= 1.54 x 10®cm“^), 
there is a single energy-momentum curve that, once the spin-polarized configuration emerges for ng > n™®®*, splits 
into three curves, one of them above and another one below the unpolarized dispersion relation. In Fig. HKb) the 
concentration is just above the threshold, so that both polarized concentrations exist below their Fermi level, which 
sits slightly higher than its unpolarized counterpart. In Fig. HKc), the curve-splitting is more significant. The topmost 
dispersion relation, corresponding to spin-up electrons, is almost entirely above the spin-polarized Fermi level, which, 
in turn, is much higher than the unpolarized level. In both Fig. HKb) andlDJc), the spin polarization tapers off at high 
wavevectors/energies, indicating that the effect weakens with carrier energy and concentration. 
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The emergence of density-dependent spin-polarized configurations is in agreement with previous experimental ob¬ 
servations of a small conductance plateau around O.SGo = jh that appears in quantum wires for specific electron 
concentrationsi^ril. According to our model, this plateau arises from the gap between the energies E {k, a) of spin- 
up and spin-down electrons in the polarized configurations, as shown in Fig. 3) Indeed, when the concentration is 
increased past the polarization threshold, e.g. by changing a gate voltage in a GaAs/AlGaAs heterostructure, one of 
the spin channels will gradually close, leading to a narrow step at e^/h instead of 2e^//i. The predicted spin-polarized 
regime is an excited state, though, and the system is expected to fall back to the unpolarized configuration (which has 
a lower average energy; see Fig. EJb)) as the density is increased, especially at low temperatures. The fact that the 
observed jh feature is narrower than the 2e^//i plateau and is enhanced with increasing temperature is consistent 
with previous conductance measurement a^^d^ . 

The theoretical results also support the interpretation of a spin-polarized excited state as the origin of the anomalous 
conductance structure at 0.7 GO observed in quantum point contacts (QPGs}^d^d^. Since the spin-polarized regimes 
occur in the high-energy excited states, the QPG potential barrier for these states is higher than that of the unpolarized 
configuration, reducing the conductance for those spin-polarized electrons. Furthermore, the electron concentration 
in spin-polarized configurations can be lowered by changing the gate voltage in the QPG heterostructure, inducing a 
further reduction of the conductance. These two facts are consistent with the temperature-enhanced pinning of the 
conductance near O.TGq at specific gate voltages. 




Figure 5. Confinement dependence of the spin-polarized regimes. Concentration no at the onset of polarization (solid 
line) and when full polarization is achieved (dashed line) as a function of: (a) hjOy for the case hjOy = hwz = hw, (b) ujy/ujz, 
with huoz = 2meV. Bx = 0 and T = 0 in both cases. Insets: close-ups of regions hio < 0.3 meV (in (a)) and ojyfwz < 0.25 (in 
(b)) showing the maximum concentration nc for which only one subband is occupied (dash-dotted line). 

As mentioned before, the concentrations at the onset of polarization (n™®®*) and at full polarization (uq”^*) also depend 
on the strength of the transverse confinement along the y- and z-directions, i.e. on the size of the cross-section of the 
wire. This is entirely due to the exchange interaction, since Eq. ([2]) predicts a confinement-independent threshold 
onset Jqj. non-interacting electrons (Cab (p) —t 0). Figure 31(a) displays the values of n™®®* (solid line) and 
(dashed line) as a function of the confinement strength hwy for a wire with a circular cross section [hiOy = hujz = huj) 
and at = 0. If fia; increases, both and increase, although the rate of increase slows down and becomes 

almost constant when Huj > 2meV. Furthermore, the range of no for which partial polarization exists (0 < |P| < 1, 
corresponding to region (III) of Fig. UKa)) becomes wider when the confinement is stronger. These data indicate that 
the wider the wire, i.e. the weaker the confinement, the lower the concentration for which polarization is possible, 
provided that the system remains in the extreme quantum limit. To satisfy this condition, the concentration must be 
smaller than a critical value nc = •y/2m* (Hujy), since hujy is the energy above the bottom of the first subband at 

which the second subband will start to become populated (when ojy < ujz)- The inset shows a close-up of the n™®®* 
and curves for small values of hut, as well as a plot of hq (dash-dotted line). For huj < 0.02meV, n™®®* > nc 
and the system is outside the scope of the model. 

The effects of asymmetric confinement on rig"®®* and rig"^* are portrayed on Figure[5l(b) as a function of the confinement 
strength ratio ujyjujz for huiz = 2meV, in the absence of a magnetic field. Both rig"®®* and rig"^* increase with the 
ratio, very rapidly for uiy < Oz and at a much slower rate for uiy > Uz- This reflects the fact that the overlap function 
Cab grows with stronger confinement (see Supplementary Figure ISTI) . For instance, for huy = Huiz = lOmeV, n™®®* = 
2.34 X 10^ cm“*, whereas for hiOy = lOmeV and hWz = 2meV (corresponding to ojyjwz = 5), rig"®®* = 1.87 x 10® cm“*. 
As in Figure [nja), the spin-polarized configurations are accessible until n™®®* > nc; as shown in the inset of Figure 
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lib), this happens when ujyluJz < 0.068 or hujy < 0.14 meV. 

Our model shows that electron-electron interactions in quasi-lD systems with longitudinal magnetic fields lead to 
complex interweaved spin-polarized regimes which differ from previous works^i, e.g. magnetic configurations oppo¬ 
site to those expected from pure Zeeman splitting in low-concentration quantum wires. These new results are the 
consequence of considering an exact 3D 1/r-Coulomb potential as opposed to a Dirac point potential^ii^, which 
enables us to include the interaction between electrons with parallel spins (prohibited with a Dirac interactioii^i) 
in the energy-momentum relation, Eq. (HD. Specifically, at zero B-field, our theory shows that two symmetric and 
degenerate spin-polarized configurations exist for all values of carrier concentration (and Fermi energy) above a con¬ 
centration threshold, even in the absence of spin-orbit interaction. The existence of the spin-polarized excited state 
in quasi-lD systems, more readily accessible as temperature is increased, is consistent with the emergence of a narrow 
conductance step near 0.5 x 2e^//i at low temperatures (^ 50 — 100mK)i2r— , and its subsequent evolution into a 
broader feature near 0.7 x 2e^//i as the temperature is increased to a few kelvins in both quantum wires and quantum 
point contact3^ d3d4,32 ^ collective spin states described in our model also account for the observations of separate 
spin modes in coupled quantum wires that were originally attributed to Luttinger-liquid behavio r^^d^ . Our findings, 
however, result from the nonlinear nature of the energy dispersion relation as the particle concentration increases, 
thereby outlining the limitations of the low-energy, linear-dispersion Luttinger model^. Such limitations are revealed 
through the observation of fast energy relaxation of particles in quantum wires, a process forbidden by Luttinger-liquid 
theory^. Moreover, our results do not conform to the Wigner-crystal picture^i^ either, since it precludes the presence 
of spin-polarized states in the absence of a magnetic field. Indeed, in this picture, the condition for polarization is 
that the Zeeman energy \Uz\ = gfiBB/2 exceed the exchange energy |f/exch|^, in contrast to our findings. 

Furthermore, our theory highlights the importance of a full 3D approach to account for the sensitivity of the spin- 
polarized configurations to varying confinement strength and asymmetry. In this context, our model uses a pair of 
parabolic potentials for transverse confinement, which is suitable for elliptical wire cross-sections. It is, however, 
general and valid for other geometries, for which the overlap function Cabip), critical in setting the concentration 
threshold Uq"®®* for the different spin-polarized regimes (Eq. (|3|)), would need to be evaluated with the corresponding 
wave functions. This is particularly relevant for spintronics applications, as it enables the design of quantum wires 
with specific spin-polarization characteristics by changing the shape and confinement of the ID channel. For instance, 
in a GaAs/AlGaAs heterostructure, the desired range of polarizations can be set by choosing a suitable physical 
separation between the split gates and a specific acceptor density in GaAs. Subsequent fine-tuning of the polarization 
could be achieved by simply adjusting the split-gate potential bias. 

Interestingly, our prediction that spin polarization is possible in symmetric quantum wires contrasts with the findings 
of Debray et al2^, whose experiments observed a spin-polarized current in a QPC only when the confinement potential 
was strongly asymmetric. It should, however, be pointed out that their QPG device is made with InAs, a material 
characterized by a strong intrinsic spin-orbit interaction. For this reason, the authors of that work attribute the 
emergence of spin polarization to a lateral spin-orbit coupling (which is a function of both the electron momentum 
and the confinement potential), as opposed to the Rashba spin-orbit interaction (which only depends on momentum). 
The extension of our model to incorporate both types of spin-orbit couplings is a topic of future investigation. However, 
the main objective of this work was to show that, even in the absence of such couplings, a variety of magnetic phases 
is achieved in quasi-lD systems by electric manipulation alone. 


METHODS 


We calculate the ID energy-momentum relation within the unrestricted Hartree-Fock approximation, considering only 
the extreme quantum limit (i.e. when only the lowest-energy subband of the confinement potential is occupied), as 
described in the Supplementary Methods. All numerical calculations (e.g. solving for Aua- in Eq. ([2)) were done 
using MATLAB. 
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SUPPLEMENTARY METHODS 

Variational calculation of the ID energy-momentum relation 

In a longitudinal magnetic field B = B^x, with the gauge set to A = B^yz, the Schrodinger equation within the 
unrestricted Hartree-Fock approximation^! reads 




1 d 

C/conf {y, z) + -m*ujgy'^ + ihhjBy-^ + C 4 i (fO 

2 az 


V'li} (^) 

t/exch [V’{i} (^] + Uziiii} (f) = (r) (SI) 


Here, {i} = {ixAyAz, (Xi = ±^} are the quantum numbers associated with the eigenenergies ujb = ^§r is the 
cyclotron frequency, q is the electron charge, m* is the effective mass, Uconi {y, z) = ^m*ujyy^ + ^m*u!'^z'^ is the lateral 
confinement potential, Uz = gy.BBxO' is the Zeeman energy term, g is the effective electron ^-factor, g.B = is ths 
Bohr magneton in the wire, and Uei and f7exch are the Hartree and exchange terms, respectively!!. 

In the extreme quantum limit, i.e. when only the lowest-energy subband of the confinement potential is occupied, we 
take the expectation value of the left-hand side of Eq. (EH), using the trial wavefunction 


V'fe. = 


1 




y/ Lx 




.1/4' 


V 


/2 


(S2) 


where a = y VlAht a^d & = y ■ In the absence of electron-electron interactions and for zero magnetic field, this 

wavefunction corresponds to the exact ground state of Eq. EH)- Then, one obtains the following expression for the 
single-particle energy: 


E {kx 3 — 2 4; 2 2 ktlOz q ^ gg^BBxCf -t- Uq\ [uq] -t- Ifexch {kx j ^) 


(S3) 


The expectation values Uei and C/exch ar( 


,27 


C^el [no] = 


(0 )^ 

167re 


-no 


Ifexch {kxj 


g 2 r +00 

g2^2g J ^P^ - kx) frlE {Px,cr)] 


where no = ^ X] It [E {p, cr)] = ^ no- is the total electron density in the wire; /y [E] = {I -|- exp [{E — p) /fcsT]} ^ 

<7 CT 

is the Fermi-Dirac distribution for a chemical potential /r; and C,ab (p), a form function which is specific to the wave- 
function shown in Eq. (IS2p . is 


, , , „/ON , <exp(-t^/2) 

Ca6(p)=8(-j/ dt -^= 


'E 


(S4) 


Cab ip) is a monotonically-decreasing function of p and it increases with increasing a and b (i.e. with stronger confine¬ 
ment strengths). A plot of Cab (p) for different values of a and b is shown in Supplementary Figure !^ 

At zero temperature, p is equal to the Fermi energy Ef and fr (E) = 9 (Ef — E (k)) = 9 {kfy^y — |A:|), where fc/(a) is 
the spin-dependent Fermi wavevector such that E (fc/(o-)) = Ef. Then, the exchange term turns into 

g2 r+kB<7)+kx 

Uexch {kx: ^) ~ ~ I dp t^ab ip) 

oZTT e 
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Calculation of Aua-, Wq"®®* and 


At T = 0, when the concentrations of both spin-up and spin-down electrons are positive, they are directly proportional 
to a spin-dependent Fermi wavevector via the relation = -Kn^. Thus, the Fermi energy Ej = E ikf(a)) be 
written in terms of Ua and np: 


Ef (n^) = 


—— + -nujy + -tu^, + -l — 




q'^Cab (0 ). 

IOttc 


-no 


327r2e , 


huJB 

r+Tmo/2+Trna- 

—Tmo/2+Trna- 


dp Cab (p) 


Ufj can be written in terms of the average concentration no/2 and the deviation from this average value = Uo-—no/2 
(so An_cr = —AtIo-). Then, since Ef is spin-independent, we get Ef {jia) = Ef {ji-a)- From this equation one obtains 
the following expression, which corresponds to Equation @ in the main text: 


TT^ * QqB. 

—noAua + a 


q^m* 


2h 


L 


27r(no/2+Ano-) 


27 r(no/ 2 —Ario-) 


dp Cab [jp) = F (An^.; no) = 0 


(S5) 


For non-interacting electrons (i.e. when the overlap function Cab {p) vanishes), the magnetic field induces a Zeeman 
splitting between spin-up and spin-down electrons. For uq > ub = tbe (unique) solution to Eq. (jS5p reduces 


to n^ = ^ 


1 - 


(^) sign(5B^CT) 


. However, for no < ns, if sign (gB^) > (<) 0, this last equation yields n.|-^) < 0, 


which must then be set to zero, so the wire exhibits full spin-down (up) polarization (n^^d) = no). 

In the general case of interacting electrons {Cab (p) 7 ^ 0), and if the concentration is low, Eq. (|S5|l yields only one 
solution for Ano- in terms of uq. However, if uq exceeds a certain threshold value ng"®®* then two additional solutions, 
or spin configurations, are possible. This is illustrated in Supplementary Figure IS^ which is a plot of E {Ana] no) 
(the left-hand side of Eq. (IS5I1 1 versus An-f- for a few values of no. As in Eigure[T]in the main text, = IT and 
hbJy = hbJz = 2meV, yielding n™®®* = 1.94 x 10®cm“^. The solutions to Eq. (ISSl) correspond to the crossings of the 
curves with the x-axis (dotted line). When the concentration is below n™®®* (dashed line), where is a single positive 
solution for An-f-, which corresponds to the “t” configuration in Fig. [TJ Meanwhile, when no > n™®®* (dash-dotted 
line), two additional negative solutions are possible; they correspond to the and configurations of Fig. [T] 
For no = n™®®* (solid line), the slope of the function F (An-f-; no) is zero at the value of An-;- for which the negative 
solution emerges. Thus, n^®®* and the corresponding An^ are the solutions of the pair of equations E (An^; no) = 0 


and 


5F(Ano-;no) 

OArio- 


= 0 . 


In the limit |Ano-| no, one can expand the integral in Eq. (1851) as a series in Ano-. Up to the third order, this 
results in a cubic equation: 


^ qgBa 

+ a— 


q^m* 

16//W 


nAnaCab (Trno) + ^ (^rAn^)^ Cab i'^Bo) 


= 0 


Here, Cat i^^o) = 


d^Cab(p) 

dp‘ 


J P—TTUq 

Ba; = 0 (Eq. ([3]) in the main text): 


This leads to the approximate expression for the non-zero values of An^. when 


An„ 


±-1 


Cab i'^Bo) 


a*no — 


Ca b (tTUq) 
27r^ 


where a* = is the effective Bohr radius. However, this solution is only valid if the argument inside the square 

root is non-negative, which means that the minimum concentration ?^o(^®*Lo) emergence of the non-zero solution 

must satisfy the identity 


Cab (^no7£Lo)) 


''0(Sx=0) 


27ra* 
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As no increases above np"®®*, An^ will continue to increase until the wire is completely spin-polarized at a concen¬ 
tration Setting Aria = ±^^ 0/2 in Eq. (IS5l) gives the following equation for the value of no for which full 

spin polarization is achieved: 


TT 


2 


4 



2 * 


-sign (crAn,,) —— 



Figure SI. (^ab (jp) as a function of p for different values of a and b. 



Figure S2. F (An-|-; no) vs. An^ for different no, with = 1 T and hwy = fiwz = 2meV. 












